NASA TECHNICAL NOTE 


CM 

vO 


i 



on 


<S 



lOAN a C F °Wu\d|gU N 

K1RTL.AND AFB, N. 


a 

b - 1 
LU 
1 1 1 

It* 

nj 

□ 



DYNAMIC RECONSTRUCTION ERRORS 
IN DIGITAL-TO-ANALOG SYSTEMS 
WITH BIOMEDICAL APPLICATIONS 

by William P. Dotson , Jr. 

Manned Spacecraft Center 
Houston , Texas 77058 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION * WASHINGTON, D. C. * APRIL 1971 


•A'* ' , ■. ■ 

, v^ V ^ \? 

^ '■ 




■ 


■»? 



TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



0133120 


1. Report No. 2. Government Accessior» No. 

NASA TN D -6296 

3. Recipient's Catalog No. 

4. Title and Subtitle 

DYNAMIC RECONSTRUCTION ERRORS IN DIGITAL-TO-ANALOG 
SYSTEMS WITH BIOMEDICAL APPLICATIONS 

5. Report Date 

April 1971 

6. Performing Organization Code 

7. Author(s) 

William P. Dotson, Jr. , MSC 

8. Performing Organization Report No. 

MSC S-254 

10. Work Unit No. 

039-00-00-00-72 

9. Performing Organization Name and Address 

Manned Spacecraft Center 
Houston, Texas 77058 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 

16. Abstract 


Various digital-to-analog conversion techniques are analyzed in order to determine their instan- 
taneous error performance. The analytical technique is curve fitting of polynomials of zero to 
third order to the given digital sample points. Experimental results are also given on digitized 
electrocardiogram data processed through a third-order polynomial curve-fitting routine. This 
report demonstrates that a trade-off of system complexity (order of the polynomial) and trans- 
mission bandwidth is possible, with no decrease in the quality of the reconstructed data. 


17. Key Words (Suggested by Author(s) I 18. Distribution Statement 

' Digital-to -Analog Converters 

' Digital-to- Analog Conversion Errors Unclassified — Unlimited 

‘ Minimum Bandwidth Systems 
‘ Interpolation 

' Telemetry Systems 


19. Security Ciassif. (of this report) 

20. Security Ciassif. (of this page) 

21. No. of Pages 

22. Price* 

None 

None 

43 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 



DYNAMIC RECONSTRUCTION ERRORS IN DIGITAL-TO- ANALOG 

SYSTEMS WITH BIOMEDICAL APPLICATIONS 

By William P. Dotson, Jr. 

Manned Spacecraft Center 

SUMMARY 


Available documentation on telemetry systems used to handle analog parameters 
from the Apollo spacecraft indicates that digital sampling rates of approximately 5f c 

to 10f c are used. (The term f £ is considered to be the highest frequency of interest 

in the spectrum of the analog signal. ) The digital samples are then transmitted to a 
sample- and- hold digital-to-analog converter where the analog signal is reconstructed. 
This report presents an analysis of the worst-case dynamic errors of this reconstruc- 
tion process and also the reconstruction processes up through the third order. The 
conclusion of this paper — based on experimental, as well as theoretical, results — 
is that certain classes of analog signals (for example, electrocardiogram signals) may 
have their sample rates reduced by a factor of 2 or more, while simultaneously the ac- 
curacy of the reconstructed signal is increased. This increase in accuracy is accom- 
plished by using higher order reconstruction processes. 

Reducing the sample rate of a parameter by a certain factor reduces the band- 
width requirement for that parameter by the same factor. Such a reduction in the 
bandwidth requirement would result in a significant easing of the present bandwidth 
problems between the remote sites and the Mission Control Center at the NASA Manned 
Spacecraft Center. Such benefits of a reduction in the sample rate (for example, eas- 
ing bandwidth problems) are obvious. Other benefits would result from increasing the 
power spectral density of a radio-frequency down link (by decreasing the bandwidth), 
without increasing the spacecraft's transmitter power. This increase in the power 
spectral density should prove of particular benefit for lunar and deep- space probes 
where the down-link signal margins are seriously low. In situations where bandwidth 
is not a serious constraint, the sample rate could remain at its present level, and the 
accuracy of the reconstructed data could be significantly improved by going to high- 
order reconstruction processes. 


INTRODUCTION 


This study arose in connection with a problem concerning electrocardiograms. 
Briefly, the problem involved satisfying the requirement to transmit an electrocardio- 
gram (EKG) for each of three astronauts from the spacecraft to the Mission Control 



Center at the NASA Manned Spacecraft Center in real time. The bandwidth necessary 
to fulfill this requirement is the product of three electrocardiograms, 200 samples/ 
sec/EKG, and 8 bits /sample, which gives 4800 bits/sec. 

Additional bits must be added to this requirement of 4800 bits /sec for overhead 
in the telemetry formats. Furthermore, it would be most desirable to be able to 
satisfy other telemetry requirements at the same time, thereby adding still more to 
the required bandwidth. 

The essence of the problem is that a bandwidth considerably in excess of 
4800 bits/sec is needed to satisfy the bandwidth requirement if the present system 
concepts are used. Unfortunately, however, the available bandwidth from a remote 
site to the Goddard Space Flight Center is limited to 4800 bits/sec, including overhead 
in the telemetry formats. 

It is clear, then, that the present system concepts must be changed to satisfy 
the basic requirement of transmitting three electrocardiograms in real time. Two 
approaches are immediately obvious: 

1. Truncate each sample from an 8-bit to a 6-bit representation. This trunca- 
tion will reduce the EKG bandwidth requirement to 3600 bits/sec at the expense of an 
increase in quantization noise. 

2. Use a redundancy removal scheme (such as a zero-order predictor). Use of 
such a scheme may or may not reduce the bandwidth requirement, depending on the 
activity of the data; and implementation of such a scheme is fairly complex. 

A third possible approach to the situation exists. If the sample rate is cut in 
half, the EKG bandwidth requirement drops to 2400 bits/sec. This sample- rate reduc- 
tion would be simple to implement (by omitting every other sample), but would result 
in a decrease in the accuracy of the reconstructed signal if a zero-order reconstruction 
process were used. However, analysis of the error performance of different- order re- 
construction processes shows that the error performance of a third-order digital-to- 
analog converter (DAC) at 100 samples/sec is better than that of a zero-order DAC at 
200 samples/sec for any spectrum component less than 20 hertz. Most of the energy 
in an EKG spectrum lies below 20 hertz, and the fastest part of the analog waveform 
(the R- segment) exhibits a fundamental frequency of approximately 17 hertz. 

The author would like to acknowledge the assistance of F. J. Wancho, who wrote 
the computer programs used in this study. 


SYMBOLS 


A amplitude 

B, C, D, E coefficients of polynomial terms 
f frequency, Hz 
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cut-off frequency, Hz 


f 

c 

f sampling frequency, samples/sec 

s 

| G(f ) | root- mean- square spectrum of g(t) 

g(t) time function 

m ratio of the spectrum cut-off rate of the signal to 6 dB/octave 

p normalized sample rate, samples/cycle (or f /f) 

T period, 1/f, sec 

t time, sec 

y(t) time function 

e error, percent or pulse- code- modulation (PCM) counts 

instantaneous error measured near the base-line crossing, percent 

^ instantaneous error measured at the peaks, percent 

t sample period, sec 

u) radian frequency, 2nf, rad/sec 

NYQUIST SAMPLING THEOREM 


The Nyquist sampling theorem states that if the root- mean- square spectrum 
| G(f)| of the time function g(t) is identically zero at all frequencies greater than f , 

then g(t) is uniquely determined if its ordinates at a series of points spaced 
l/(2f c ) seconds apart are known, with the series extending throughout the time do- 
main (ref. 1). The theorem, as stated, does not consider two problems that always 
occur in the real world. The two problems are related to the following two statements. 

1. Ideal band-limited signals do not exist. 

2. Even if an ideal band-limited signal did exist, it would be necessary to pass 
its sampled values through an ideal low-pass filter with a cut-off frequency of f /2 in 
order to reconstruct exactly the analog signal. 
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Downing (ref. 1) has shown that the real-world signal may be considered to be 
band limited at the negative 3- decibel frequency for system design purposes, provided 
the following factors are known: 

1. The rate of spectrum cut-off of the analog signal 

2. The acceptable root- mean- square aliasing error 

Knowledge of these factors permits a 

choice of the minimum sample rate to the 

negative 3- decibel frequency ratio. In the 

curves shown in figure 1 (ref. 1), m is 

the ratio of the spectrum cut-off rate of 

the signal to 6 dB/octave. It should be 

noted in using the curves shown in figure 1 

that they presuppose the use of an ideal 

low-pass filter with a cut-off frequency of 

f /2 for the digital- to- analog conversion 
s 

process. It should also be noted that the 
curves in figure 1 give the root- mean- 
square error and not the instantaneous 
error at any particular point on the ana- 
log waveform. The instantaneous error 
may be of considerable importance in the 
analysis of particular classes of data and 
is, in general, much greater than the 
root- mean- square error. 

Ideal or optimum filtering for digital-to-analog conversion (or high-order recon- 
struction processes) minimizes both the root- mean- square and the instantaneous 
errors. Some remarks on optimum filtering for digital-to-analog conversion may be 
found in reference 1. This report is concerned with the instantaneous errors of digital- 
to-analog conversion processes by which a curve described by a polynomial in t is 
fitted to the given sample values. 



Normalized sampling frequency, f s /f c 

Figure 1. - Root- mean- square aliasing 
error as a function of normalized 
sampling frequency. 


DYNAMIC-RECONSTRUCTION ERRORS 


The instantaneous error as a function of time between any two arbitrarily placed 
sample points is both difficult to find and difficult to use. Accordingly, in this section, 
the instantaneous errors will be found in two regions of interest. These regions are 
(1) near the base-line crossing of the analog signal and (2) at the peak of the analog sig- 
nal. The time ordinates of the given samples have been selected so that both a mean- 
ingful result and a straightforward analysis can be obtained. Different choices in the 
placement of the sample points will, of course, result in different error equations, 
which, in most instances, will result in lower errors. The implication is that the 
analysis given is of a worst- case nature; and if random starting points for the sample 
locations were assumed, an average error, over many cycles, of approximately half 
the value of the worst- case errors could be expected. 
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Zero-Order DAC 


Figure 2 is used to find the instanta- 
neous error midway between two sample 
points for a zero- order DAC. This error 
is near the base-line crossing. The in- 
stantaneous error at t = r/2 for the 
sample-point locations shown in figure 2 
is 




Figure 2. - Geometry for calculating 
for a zero-order DAC. 


y(t) = A sin (cot) 


(lb) 


y' (t) = A sin [ co (0)] 0 < t < t 


(lc) 



= A sin 



(Id) 


= 0 


y = A 

J max 


(le) 

(l f) 


% 


= 100 sin 


in (coj) 


2tt 


co = 


T 

T = — 

P 


(lg) 


where p = normalized sample rate in samples/cycle 


e. = 1 00 sin — 
p 


( 2 ) 
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Figure 3 is used to find the instanta- 
neous error at the peaks for a zero- order 
DAC. The instantaneous error at t = 0 
for the sample-point locations shown in 
figure 3 is 


e 

P 


y(Q) - y'(Q) 


max 


( 100 ) 


(3a) 


y(t) = A sin 



(3b) 



Time, t 


Figure 3. - Geometry for calculating 

e for a zero- order DAC. 

P 


y' (t) = A sin 


U> 


/ T T' 
1“ 2 + 4, 


i<t< - 
2 < 2 


(3c) 


y(0) = A sin 


in |o> (o + 


= A 


(3d) 


y’ (0) = A sin [w (- J 


(3e) 


e 

P 




(3f) 


(3g) 


e = 100<1 - sin 

P I 



(3h) 


By using trigonometric identities, the equation 


e = 100(l - cos 
P \ P / 

can be obtained. Equations (2) and (4) are plotted in figure 4. 


( 4 ) 
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Error, «, percent 


I 
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Normalized sample rate, p, samples/cycle 



First-Order DAC 


T - 


Figure 5 is used to find the instanta- 
neous error midway between two sample 
points for a first-order DAC. This error 
is near the base-line crossing. The in- 
stantaneous error at t = r/2 for the 
sample-point locations shown in figure 5 
is 


e b = 



(5a) 


I I I 1 llii 

0 “ t 2t 4t 6t 8-r 10 t 12t 

* Time , t 

Figure 5. - Geometry for calculating 
for a first-order DAC. 


y(t) = A sin (o>t) 


(5b) 



y'(t) = o + 


|A sin (lot )| 


0 < t < T 


(5c) 


y(j) = A sin (<s j) 
y, (j) = \ A sin (lot) 


(5d) 

(5e) 


y = A 
J max 


(5f) 



(5g) 


eb = 100(sin^-|sin^) 


( 6 ) 
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Figure 6 is used to find the instanta- 
neous error at the peaks for a first-order 
DAC. The instantaneous error at t = 0 
for the sample-point locations shown in 
figure 6 is 


e 

P 


y(Q) - y'(Q) 

y 

•'max 


( 100 ) 


(7a) 


y(t) = A sin 



(7b) 



Figure 6. - Geometry for calculating 

e for a first-order DAC. 

P 


y' (t) = A sin 




A sin [c 

4*11 

- A sin 

[4W2 

1 



T 



t , r 
2 ~ 1 < 2 


(7c) 


y(0) = A sin (w ^ = A 
y'(0) = A sin ^ +0 


(7d) 


(7e) 


e 

P 


y = A 
J max 


lOojl - sin [a>(- 1 + J) | 



(7f) 


(7g) 


e 

P 


100 


I- 



(7h) 
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By using trigonometric identities, the equation 


e p = 100(l - cos (8) 

can be obtained. Equations (6) and (8) are plotted in figure 7. 



40 30 20 10 7 5 4 3 2 

Normalized sample rate, p, samples/cycle 

Figure 7. - Error as a function of normalized sample rate, first-order DAC. 
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Second-Order DAC 

Figure 8 is used to find the instanta- 
neous error midway between two sample 
points for a second- order DAC. This 
error is near the base-line crossing. The 
instantaneous error at t = 0 for the 
sample-point locations shown in figure 8 
is 


= y(°) - V(0) (io 0 ) (9a) 

^max I 

y(t) = A sin Ut + '-j (9b) 

y' (t) = Ct 2 + Dt + E 



^tOj. j: t. h ^ n- Ur 

22222 222 
Time, t 


Figure 8. - Geometry for calculating 
for a second-order DAC. 

- 5 * t < % ( 9c ) 


y(0) = A sin (ol> ^ (9d) 

y’(0) = E (9e) 

W = A (9f) 

e b = 100[sin(uil)-|] (10) 

To find E, solve the set 

y (- 1) = c (- 5) 2 * D (- 5) - E = A stn [" (- 5 * 5)] = 0 (lla) 

2 

y '( J) = c(l) + d(J) + E - A sin (an) (lib) 
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y '(^ =c (r) + D (^-) + E = A sin ( 2wt ) 


(He) 


E = - g- [ sin (2arr) - 6 sin (wr)] 
= 100 1 sin ^co ^ [6 sin (ut) - sin (2ojt)] 


2tt 


CO = 


T 

T - — 

P 


% 


= 100 


77 1 

sin — „ 6 sin 

p 8 


2tt . 47t\ 

sin — 

P P/ 


Figure 9 is used to find the instanta- 
neous error at the peaks for a second- 
order DAC. The instantaneous error at 
t = 0 for the sample-point locations shown 
in figure 9 is 


e = 

P 


y(o) - y'(0) 


max 


y(t) = A sin 


y’ (t) = Ct + Dt + E 


(100) (13a) 


(13b) 


co t + 


1 < t < ll 

2 < 2 


(lid) 


(lie) 


(12) 



(13c) 

y(0) = A sin (w ^ = A 
y'(o) = E 


Figure 9. - Geometry for calculating 

e for a second- order DAC. 

P 


y = A 
max 


(13d) 

(13e) 
(1 3f) 
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( 14 ) 


e = 100 
P 




To find E, solve the set 


y ' (- 5 ) = C (- if + D (-i) + E = A Sin [“ (- 5 * ?) 


Hi) = c (if - °{i ) - E = A sl " Hi * ?) 


y, (x> = c fir) 2 + D (ir) + E - A sl " [«(t + j) 


E = ^{3 sin 


T T N 
ft,l “ 2 + T 


j + 6 sin L(l + 1) 


- sin 


, 3t T 

2 + 4 


277 


OJ = 


= I 

7 P 


e p = 100 1 1 - sin 


(tt)! + 3 sin Rtf) 


- sin 


p + 6 
2p 


By using trigonometric identities, the equation 


e = 100 
P 


, 1 /_ 77 77 \ 

1 - O 9 cos — cos 
8 V p 3p/ 


can be obtained. Equations (12) and (16) are plotted in figure 10. 


(15a) 

(15b) 

(15c) 

(1 5d) 
(15e) 

(16) 
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Third-Order DAC 

Figure 11 is used to find the instan- 
taneous error midway between two sample 
points for a third- order DAC. This error 
is near the base-line crossing. The in- 
stantaneous error at t = 0 for the 
sample-point locations shown in figure 11 
is 




y(o) - y’(o) 


max 


( 100 ) 


(17a) 



y(t) = A sin 



(17b) 


Figure 11. - Geometry for calculating 
for a third-order DAC. 


y' (t) = Bt 3 + Ct 2 + Dt + E 


_ < t < 

2 ~ 1 < 2 


(17c) 


y(0) = A sin 



y’(0) = E 

y = A 

max 



To find E, solve the set 


( 3t\ 

( 3t\ 3 „/ 

3t\ 2 , 1 

' 3r\ _ . , 

(t \~ 

\ 2/ 

‘ B \- t) +c( 

nr) +D > 

r + E = A sin 



y '(~ §) = B (" i) +C ("J) + D ('s) +E ' Asln ("y)=° 


(17d) 

(17e) 

(17f) 

(18) 

(19a) 

(19b) 
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-\ 3 h 

-\ 2 h 

\ 

/T \" 

)- B (; 

i) +c (; 

0 +d G 

-J + E = A sin 

"(2 + T ) 


(19c) 


y’(r) = b (t ) 3 + c (t ) 2 + d (t) + E * A sin [“(I - 2t )] <19d) 


E 


= ^| 9 sin [o>(| + r)] - sin o>(£ - t)~ - sin «(£ + 2r)"| 


(19e) 


- sin 


i T „ > 

w( .2 +2t J 


( 




<9 sin 

"(2 + T ) 

- sin 

"(2- t ) 


u> = 


T = 


27 T 
T 

T 


(19f) 


= lOO^sin 


P + 1 


1 

16 


9 sin 


77^ 

V p 


sin 


P- 2 


- Sin 


p + 4 
P 


(19g) 


By using trigonometric identities, the equation 


e, = 100 
b 


. 77 1 /. . 277 ^ . 477 

- sin - + tt* ( 1 0 sin sin — 

P 16 V P P 


( 20 ) 


can be obtained. 

Figure 12 is used to find the instan- 
taneous error at the peaks for a third- 
order DAC. The instantaneous error at 
t = 0 for the sample-point locations shown 
in figure 12 is 


e 

P 


y(Q) - y'(Q) 


max 


( 100 ) 


(21a) 



y(t) = A sin 



(21b) 


Figure 12. - Geometry for calculating 

e for a third-order DAC. 

P 
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(21c) 


y’ (t) = Bt 3 + Ct 2 + Dt + E “ IT ~ 1 < T 


y(0) = A sin 



= A 


y'(0) = E 

y = A 
max 

£p = 100(1 . f) 


To find E, solve the set 


y* 


( 3rY 

3 n( 3t\ 2 

3tN 


/T 

3t\ 

\ 2J 

+ c (r“27 +D v 

2 J 

+ E = A sin 

L w W- 

2/ 


y’(- \)= B (- if + 4 if + D (- i)+ E “ A s,n [“(? - 0 


/, 

r\ A, 

.\3 / 

A 2 / 

A 


y’G 

44 

s) +c (l 

s) + 4 

i ) + E = A sin 

w (t + 2) 



(21 d) 

(21 e) 
(21f) 

( 22 ) 

(23a) 

(23b) 

(23c) 

(23d) 

(23e) 
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By using trigonometric identities, the equation 


e p = ioo [ i ■ K 9 cos l ■ cos t) (24) 

can be obtained. Equations (20) and (24) are plotted in figure 13. 



40 30 20 10 7 5 4 3 2 

Normalized sample rate, p, samples/cycle 

Figure 13. - Error as a function of normalized sample rate, third-order DAC. 
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EXPERIMENTAL RESULTS 


As a means of checking the validity 
of the results of the preceding section, 
several variations of the data processing 
technique were simulated on a digital PCM 
system with a hybrid computer program. 
Conceptually, the system worked as indi- 
cated in figure 14. Table I is an explana- 
tion of the stages in figure 14 and in 
tables n and ni. 


Stage A Stage B Stage C Stage D 



Figure 14. - Block diagram of the test 
system. 


The experiment- control data processing technique was chosen to be the follow- 
ing: Stage A; stage B, items 1 and 3; stage C, item 1; and stage D. This variation 
is technique 08. The variations that are presented in the form of strip-chart record- 
ings in figures 15 to 24 are given in table II. Each figure (figs. 15 to 24) shows the 
same set of data as it would appear after processing under the six variations listed in 
table II. 


Tables IV to Vin are condensed summary listings of normal sinus data; the 
tables show peak values for the P-, Q-, R-, S-, and T- segments of the analog wave- 
form on a cycle-to- cycle basis as the values would appear under the different proc- 
essing techniques. Amplitude values are given in terms of PCM counts, ranging 
from 0 to 255. The data-processing-technique variations that were compared are 
shown in table III. 


The analog data were artificially generated by a polyrhythm generator so that 
the waveform would remain stable from cycle to cycle. Hence, any variations in peak 
values of the output (the condensed summary listings) are caused by errors introduced 
by the different data processing techniques. 


Figures 25 to 29 are graphic presentations of the results found in tables IV to 
VIH. These figures show the peak errors e and the average peak errors, average 


<y for the P-, Q-, 


P 

R-, S-, and T- segments of the data as a function of the data proc- 


essing technique. Also shown, in the interest of systems- engineering requirements, 
is the bandwidth (in bits/sec) required for transmitting the same set of data (three 
electrocardiograms) for each data processing technique. 


The experimental results (figs. 25 to 29) confirm the analytical results (figs. 4 
and 13). In performing a numerical comparison of the experimental results and the 
analytical results, the following factors must be considered: 


1 . Because of the digital nature of the system, the experimental data contain 
quantization noise, whereas the analytical results were derived for analog data without 
quantization noise. 

2. Because of the analytical base (sinusoids), the error definition for the analyt- 
ical results is in terms of percent of the maximum amplitude, whereas the experimental 
results, based on complex waveforms, required an error definition in terms of percent 
of full scale (255 PCM counts). 
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TABLE I. - DESCRIPTION OF TEST-SYSTEM STAGES AND CODES FOR 
STRIP-CHART RECORDINGS AND DIGITAL DATA 
(a) Test-system stages 


Stage 

Description 

Data 

A 

Several types of EKG analog data 

1. Sinus bradycardia, 40beats/min 


were used as input for the 

2. Normal sinus, 80 beats/min 


analog-to-digital converter. 

3. Sinus tachycardia, 120 beats/min 

4. First-degree block at 80 beats/min 

5. Second-degree block at 80 beats/min 

6. Third-degree block at 80 beats/min 

7. Atrial flutter 

8 . Atrial fibrillation 

9 . Ventricular tachycardia 
10. Ventricular fibrillation 

B 

Various analog-to-digital conver- 

1. 200 samples/sec 


ter variations were used. 

2. 100 samples/sec 

3. Data quantized to 8 bits (256 levels) 

4. Data quantized to 7 bits (128 levels) 

5. Data quantized to 6 bits (64 levels) 

C 

Interpolation was performed by a 
digital computer. 

1 . If the input rate was 200 samples/ 

sec, no interpolation was per- 
formed, and the output was 
200 samples/sec (throughput). 

2. If the input rate was 100 samples/ 

sec, third-order interpolation was 
used to generate intermediate 
samples. 

a. For strip- chart recordings, an 

output rate of 200 samples/ 
sec was used. 

b. For a check against the analyt- 

ically predicted results, a 
summary listing of data 
points output at 1000 samples/ 
sec was provided by the digi- 
tal program. 

D 

The DAC was a standard zero- 
order DAC. 

— 


(b) Codes for strip-chart recordings and digital data 


Code 

Meaning 

Description 

O 

Original data at 200 samples/sec 

This was the control used in the 
experiment. 

T 

Truncated data 

The number of possible levels used to 
describe an analog amplitude was re- 
duced. This truncation was performed 
on an individual- sample basis. 

C 

Calculated data 

The implication is that only 100 samples/ 
sec of data were entered into the sys- 
tem and that these plus intermediate 
calculated samples, based on third- 
order polynomial computations, were 
output. 

6, 7, and 8 


These numbers indicate the number of 
bits used to represent the analog am- 
plitude of each sample value. 








r * 
{ . \ 

TABLE E. - 

DATA-PROCESSING- TECHNIQUE VARIATIONS 
FOR STRIP- CHART RECORDINGS 

/A , j 

; J 

Technique 

Components of the data processing technique 

If-/ . 

08 

Stage A 

' \ . ' 

^ , • ■ 


Stage B, items 1 and 3 



Stage C, item 1 

$; : 


Stage D 

; '/{{ . [ t 

T6 

Stage A 

! ' • 


Stage B, items 1 and 5 

■ v ti- *< - , 


Stage C, item 1 

% A , ' 


Stage D . 

'}y , ,, f 

C8 

Stage A 

wy; ' ' ' j * 
* ■?}{, • ‘ ' 


Stage B, items 2 and 3 



Stage C, item 2a 

; v - 4 

«, ’ -i • 


Stage D 


C7 

Stage A 



Stage B, items 2 and 4 



Stage 'C, item 2a 

' 


Stage D 

j 

a C7/T6 

Stage A 

:v ' " 


Stage B, items 2 and. 4 

, > 


Stage C, item 2a 



Stage D 


C6 

Stage A 



Stage B, items 2 and 5 



Stage C, item 2a 



Stage D 


cl 

The T6 indicates that the output data were quantized to a 6-bit level after com- 
putation and before output. 


i 

r 


21 









TABLE m. - DATA- PROCESSING- TECHNIQUE VARIATIONS 
FOR CONDENSED SUMMARY LISTINGS 


Technique 

Components of the data processing technique 

08 

Stage A, item 2 
Stage B, items 1 and 3 
Stage C, item 1 

T6 

Stage A, item 2 
Stage B, items 1 and 5 
Stage C, item 1 

C8 

Stage A, item 2 
Stage B, items 2 and 3 
Stage C, item 2b 

C7 

Stage A, item 2 
Stage B, items 2 and 4 
Stage C, item 2b 

a C7/T6 

Stage A, item 2 
Stage B, items 2 and 4 
Stage C, item 2b 

C6 

Stage A, item 2 
Stage B, items 2 and 5 
Stage C, item 2b 


3 , 

The T6 indicates that the output data were quantized to a 6-bit level after com- 
putation and before output. 
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Figure 15. - Strip-chart recordings of sinus bradycardia, 40 beats/min. 






Figure 16. - Strip-chart recordings of normal sinus, 80 beats/min. 
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Figure 17. - Strip- chart recordings of sinus tachycardia, 120 beats /min. 
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Figure 18. - Strip- chart recordings of a first- degree block, 80 beats/min 





Figure 19. - Strip- chart 
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Figure 20. - Strip-chart recordings of a third-degree block, 80 beats/min. 


28 




Figure 21. - Strip-chart recordings of atrial flutter. 
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Figure 22. - Strip-chart recordings of atrial fibrillation. 
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Figure 23. - Strip- chart recordings of ventricular tachycardia. 
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Figure 24. - Strip-chart recordings of ventricular fibrillation. 
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TABLE IV. - CONDENSED SUMMARY LISTING OF PEAK AMPLITUDES, 
P- SEGMENT, NORMAL SINUS 


Peak amplitude, PCM counts, for 
data processing technique — 



08 

T6 

C8 

C7 

C7/T6 

C6 

Cycle 1 

135 

132 

135 

134 

132 

132 

Cycle 2 

135 

132 

134 

134 

132 

132 

Cycle 3 

135 

132 

135 

134 

132 

132 

Cycle 4 

134 

132 

134 

134 

132 

132 

Cycle 5 

135 

132 

135 

134 

132 

132 

(jycle 6 

135 

132 

135 

134 

132 

132 

Cycle 7 

135 

132 

135 

134 

132 

132 

Cycle 8 

136 

136 

136 

136 

136 

136 

Cycle 9 

134 

132 

134 

134 

132 

132 

Cycle 10 

135 

132 

135 

134 

132 

132 

Average 

134. 9 

132. 4 

134.8 

134. 2 

132.4 

132. 4 

Analog 

136 

136 

136 

136 

136 

136 

Average e 

P 

1.1 

1.6 

1.2 

1.8 

1.6 

1.6 

e 

P 

2 

4 

2 

2 

4 

4 


TABLE V. - CONDENSED SUMMARY LISTING OF PEAK AMPLITUDES, 
Q-SEGMENT, NORMAL SINUS 


Measurement 


Peak amplitude, PCM counts, for 
data processing technique — 


08 

T6 

C8 

C7 

C7/T6 

C6 

Cycle 1 

117 

116 

117 

116 

116 

112 

Cycle 2 

117 

116 

116 

114 

112 

112 

Cycle 3 

117 

116 

117 

116 

116 

112 

Cycle 4 

117 

116 

116 

114 

112 

112 

Cycle 5 

116 

116 

116 

114 

112 

112 

Cycle 6 

117 

116 

117 

116 

116 

112 

Cycle 7 

117 

116 

115 

114 

112 

112 

Cycle 8 

117 

116 

115 

114 

112 

112 

Cycle 9 

117 

116 

115 

114 

112 

112 

Cycle 10 

116 

116 

114 

112 

112 

112 

Average 

116. 8 

116. 0 

115.8 

114.4 

113.2 

112. 0 

Analog 

116 

116 

116 

116 

116 

116 

Average 

.8 

.0 

.2 

1.6 

2.8 

4.0 

€ 

P 

1 

0 

2 

4 

4 

4 



TABLE VI. - CONDENSED SUMMARY LISTING OF PEAK AMPLITUDES. 


R-SEGMENT, NORMAL SINUS 


Measurement 

Peak amplitude, PCM counts, for 
data processing technique — 

08 

T6 

C8 

C7 

C7/T6 

C6 

Cycle 1 

187 

184 

186 

186 

184 

184 

Cycle 2 

187 

184 

187 

186 

184 

184 

Cycle 3 

188 

188 

188 

188 

188 

188 

Cycle 4 

187 

184 

187 

186 

184 

184 

Cycle 5 

189 

188 

189 

188 

188 

188 

Cycle 6 

187 

184 

187 

186 

184 

184 

Cycle 7 

187 

184 

187 

186 

184 

184 

Cycle 8 

188 

188 

187 

186 

184 

184 

Cycle 9 

187 

184 

186 

184 

184 

184 

Cycle 10 

189 

188 

187 

186 

184 

184 

Average 

187.6 

185.6 

187.1 

186. 2 

184.8 

184.8 

Analog 

189 

189 

189 

189 

189 

189 

Average 

1.4 

3.4 

1.9 

2.8 

4.2 

4. 2 

€ 

P 

2 

5 

3 

5 

5 

5 


TABLE VII. - CONDENSED SUMMARY LISTING OF PEAK AMPLITUDES, 
S-SEGMENT, NORMAL SINUS 


Measurement 

Peak amplitude, PCM counts, for 
data processing technique — 

08 

T6 

C8 

C7 

C7/T6 

C6 

Cycle 1 

111 

108 

111 

110 

108 

108 

Cycle 2 

110 

108 

109 

108 

108 

104 

Cycle 3 

112 

112 

111 

110 

108 

108 

Cycle 4 

110 

108 

110 

110 

108 

108 

Cycle 5 

111 

108 

110 

110 

108 

108 

Cycle 6 

111 

108 

111 

110 

108 

108 

Cycle 7 

110 

108 

110 

108 

108 

104 

Cycle 8 

111 

108 

111 

110 

108 

108 

Cycle 9 

110 

108 

109 

108 

108 

104 

Cycle 10 

111 

108 

110 

108 

108 

104 

Average 

110. 7 

108.4 

110.2 

109.2 

108.0 

106. 4 

Analog 

11 0 

110 

110 

110 

110 

110 

Average e 

. 7 

1.6 

.2 

.8 

2.0 

3.6 

€ 

P 

2 

2 

1 

2 

2 

6 
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TABLE Vni. - CONDENSED SUMMARY LISTING OF PEAK AMPLITUDES, 


T-SEGMENT, NORMAL SINUS 


Measurement 


Peak amplitude, PCM counts, for 
data processing technique — 


08 

T6 

C8 

C7 

C7/T6 

C6 

Cycle 1 

154 

152 

154 

154 

152 

152 

Cycle 2 

155 

152 

155 

154 

152 

152 

Cycle 3 

154 

152 

154 

154 

152 

152 

Cycle 4 

155 

152 

155 

154 

152 

152 

Cycle 5 

155 

152 

155 

154 

152 

152 

Cycle 6 

154 

152 

154 

154 

152 

152 

Cycle 7 

155 

152 

155 

154 

152 

152 

Cycle 8 

154 

152 

154 

154 

152 

152 

Cycle 9 

155 

152 

155 

154 

152 

152 

Cycle 10 

155 

152 

155 

154 

152 

152 

Average 

154. 6 

152. 0 

154.6 

154. 0 

152. 0 

152. 0 

Analog 

155 

155 

155 

155 

155 

155 

Average 

. 4 

3. 0 

. 4 

1.0 

3. 0 

3. 0 

e 

P 

1 

3 

1 

1 

3 

3 
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Data processing technique 


Figure 27. - Bandwidth requirement and average e p as functions of data 
processing technique for the R- segment of a normal sinus wave. 
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CONCLUSIONS 


The data processing technique using calculated data quantized to an 8-bit level 
will provide (within the limits of quantization error) as accurate a reproduction of the 
electrocardiograms as the technique using the original data quantized to an 8-bit level; 
and the technique using the calculated 8 -bit data will provide this reproduction with 
only half the bandwidth required by the technique using the original 8 -bit data. Fur- 
thermore, the technique using the calculated 8 -bit data is preferable, from both user 
and system viewpoints, to the technique using truncated data quantized to a 6 -bit level, 
which is presently being implemented; that is, the technique using calculated 8 -bit 
data is better in error performance, contains less quantization noise, and requires 
less bandwidth than the technique using truncated 6 -bit data. 

However, a problem exists in implementing the technique using calculated 8 -bit 
data. The problem is a result of the system philosophy of designing only high-speed 
data-transfer formats containing whole-number multiples of 2400 bits. The technique 
using calculated 8 -bit data requires precisely 2400 bits/sec for three electrocardio- 
grams. Thus, because of this system philosophy, no allowances are made for imped- 
ance pneumograms or for overhead (both of which are necessary) in the same format. 

The technique using calculated data quantized to the 7-bit level must then be con- 
sidered. This technique is a solution to the problem mentioned in the previous para- 
graph and, from both user and system viewpoints, is also preferable to the technique 
using truncated 6 -bit data. In an effort to reduce the required bit -rate bandwidth fur- 
ther, the techniques using calculated 7-bit data truncated to 6 bits and using calculated 
6 -bit data were tried. With these techniques, the error performance became slightly 
worse, but these techniques may be usable in the event of severe system constraints. 

Finally, it should be noted that while the interpolation in this study was done by 
a digital computer, the interpolation could also be accomplished with a small, spe- 
cialized analog computer. Such a device would generate, between the sample points, 
a piecewise continuous waveform which would effectively insert an infinite number of 
calculated points between the given sample points. Such an output would be more 
pleasing to the eye and would allow extremely accurate computation of segment dura- 
tions of the electrocardiogram waveform, since the third- order digital- to- analog con- 
verter performance is even better near base-line crossings. 
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